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Abstract 

The optimal power-flow problem (OPF) has always played a key role in the 
planning and operation of power systems. Due to the non-linear nature of the 
AC power-flow equations, the OPF problem is known to be non-convex, there¬ 
fore hard to solve. During the last few years several methods for solving the 
OPF have been proposed. The majority of them rely on approximations, often 
applied to the network model, aiming at making OPF convex and yielding inex¬ 
act solutions. Others, kept the non-convex nature of the OPF with consequent 
increase of the computational complexity, inadequateness for real time control 
applications and sub-optimality of the identified solution. Recently, Farivar and 
Low proposed a method that is claimed to be exact for the case of radial distri¬ 
bution systems under specific assumptions, despite no apparent approximations. 
In our work, we show that it is, in fact, not exact. On one hand, there is a mis¬ 
interpretation of the physical network model related to the ampacity constraint 
of the lines’ current flows. On the other hand, the proof of the exactness of the 
proposed relaxation requires unrealistic assumptions and, in particular, (i) full 
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controllability of loads and generation in the network and (ii) no upper-bound 
on the controllable loads. We also show that the extension of this approach to 
account for exact line models might provide physically infeasible solutions. In 
addition to the aforementioned convexification method, recently several contri¬ 
butions have proposed OPF algorithms that rely on the use of the alternating- 
direction method of multipliers (ADMM). However, as we show in this work, 
there are cases for which the ADMM-based solution of the non-relaxed OPF 
problem fails to converge. To overcome the aforementioned limitations, we pro¬ 
pose a specific algorithm for the solution of a non-approximated, non-convex 
OPF problem in radial distribution systems. In view of the complexity of the 
contribution, this work is divided in two parts. In this first part, we specifically 
discuss the limitations of both BFM and ADMM to solve the OPF problem. 
Keywords: OPF, ADMM, decomposition methods, method of multipliers, 
convex relaxation, active distribution networks. 


1. Introduction 

The category of optimal power-flow problems (OPFs) represents the main set 
of problems for the optimal operation of power systems. The first formulation 
of an OPF problem appeared in the early 1960s and has been well-defined ever 
since [1]. It consists in determining the operating point of controllable resources 
in an electric network in order to satisfy a specific network objective subject 
to a wide range of constraints. Typical controllable resources considered in the 
literature are generators, storage systems, on-load tap changers (OLTC), flexible 
AC transmission systems (FACTS) and loads (e.g., [2, 3, 4, 5, 6]). The network 
objective is usually the minimization of losses or generation costs, and typical 
constraints include power-flow equations, capability curves of the controllable 
resources, as well as operational limits on the line power-flows and node voltages 

(e-g., [7]). 

The OPF problem is known to be non-convex, thus difficult to solve effi¬ 
ciently (e.g., [8, 9, 10]). Since the problem was first formulated, several tech- 
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niques have been used for its solution. Among others, non-linear and quadratic 
programming techniques, Newton-based methods, interior point methods in 
the earlier years, as well as heuristic approaches based on genetic algorithms, 
evolutionary programming, and particle-swarm optimization in recent years 
(e.g., [11, 12, 13, 14]). These techniques, even though they have been shown 
to sucessfully solve instances of the non-convex OPF problem, seek to find a 
local optimal solution of the OPF. They, generally, utilize powerful general 
purpose solvers or in-house developed software but they cannot guarantee the 
identification of the global optimal solution. In general, they are character¬ 
ized by high computational complexity. The first category of approaches make 
use of gradient-based optimization algorithms or even require the use of hes¬ 
sian matrices related to the problem. Therefore, such techniques require several 
assumptions on the OPF problem formulation such as analytic and smooth ob¬ 
jective functions. Heuristics have been applied widely in the literature as a 
solution technique, for instance in cases where the OPF problem is non-smooth, 
non-differentiable and highly non-linear. 

Recently, the OPF problem is becoming more compelling due to the increas¬ 
ing penetration of embedded generation in distribution networks, essentially 
composed by renewable resources 1 . The distributed nature of such resources, 
as well as their large number and potential stochasticity increase significantly 
the complexity and the size of the OPF problem and bring about the need for 
distributed solutions. In this direction, several algorithms have been proposed 
in the literature to handle large-scale OPF problems (e.g., [15, 16, 17]). Addi¬ 
tionally, several contributions have proposed specific distributed algorithms for 
the solution of the OPF problem. In [18, 19] the authors design a dual-ascent 


Ht is worth noting that transmission and distribution systems are different with respect to 
(i) topology, (ii) electrical line parameters, (iii) power flow values, (iv) nature and number of 
controllable devices. Therefore, these systems require dedicated OPF algorithms that account 
for their specific characteristics. The focus of this work is on OPF algorithms specifically 
designed for the case of distribution networks. 
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algorithm for optimal reactive power-flow with power and voltage constraints. 
In [20, 21] dual decomposition is used as the basis for the distributed solution 
of the OPF problem. Finally, a significant number of contributions propose 
distributed formulations of the OPF problem that are based on the alternating 
direction method of multipliers (ADMM) (e.g., [22, 20, 23, 24, 25, 26]). 

Recently a lot of emphasis is put on the convexification of the OPF prob¬ 
lem. The reason behind this emerging trend is that convex problems provide 
convergence guarantees to an optimal solution and therefore such methods can 
be deployed within the context of control applications for power systems and 
specifically distribution networks. However, most of the proposed convexifica¬ 
tion schemes either do not guarantee to yield an optimal solution or they are 
based on approximations that convexity the problem in order to guarantee con¬ 
vergence. These approximations, often, either lead to (i) misinterpretation of 
the system model [27] or (ii) solutions that, even though mathematically sound, 
might be far away from the real optimal solution, thus having little meaning for 
the grid operation [28]. 

Recently, Farivar and Low proposed in [29, 30] a convexification of the prob¬ 
lem that is claimed to be exact for radial networks. In Part I of this paper, we 
show that this claim is not exact, as the convexification of the problem leads to 
an inexact system model. We also show that the method of ADMM-based de¬ 
composition, which comes together with the convexification, does not work for a 
correct system model. In this first part of the paper we focus on the Farivar-Low 
convexification and ADMM algorithms since they are considered as the most 
prominent ones by the recent literature on the subject. As an alternative, we 
propose in Part II an algorithm for the solution of the correct AC OPF problem 
in radial networks. Like ADMM, it uses an augmented Lagrangian, but unlike 
ADMM, it uses primal decomposition [31] and does not require that the problem 
be convex. We consider a direct-sequence representation of the electric distri¬ 
bution grid and we present both a centralized and a decentralized asynchronous 
version of the algorithm. 

The structure of this first part is as follows. In Section 2 we present the 
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generic formulation of the OPF problem in radial distribution systems and we 
classify several OPF algorithms based on the approximations and assumptions 
on which they rely. In Section 3 we discuss the limitations and applicability of 
the Farivar-Low formulation of the OPF problem proposed in [29, 30]. We pro¬ 
vide, in Section 4, the ADMM-based solution of the original non-approximated 
OPF problem. In the same section, we highlight specific cases where the 
ADMM-based algorithm fails to converge. Finally, we provide the main ob¬ 
servations and concluding remarks for this part in Section 4.2. 


2. Generic Formulation of the OPF Problem 

2.1. Notation and Network Representation 

In the rest of the paper, we consider a balanced radial network composed 
of buses (B). lines (£), generators (£/) and loads (C). The network admittance 
matrix is denoted by Y. Several generators/loads can be connected to a bus b€B. 
We denote that a generator g£Q or a load cGC is connected to a bus by “gGb” and 
“c€&”. We assume that the nodal-power injections are voltage-independent. A 
line £gC is represented using its exact 7r-equivalent model and it has a receiving 
and a sending end denoted by £ + and £~. Each line is connected to two adjacent 
buses: (3(£ + ) and /3(£~), respectively. Y) denotes the longitudinal admittance 
of a line, Y e + (Y^~) is the shunt capacitance at the receiving (sending) end of 
the line 2 . The notation adopted is shown in detail in Fig. 1 where the network 
branch connecting the generic network nodes i and j is represented. 


2.2. Generic OPF Formulation 

The traditional formulation of the OPF problem consists in minimizing a 
specific network objective: 


min 

S g ,S c ,S+,S^J+,If,V b 


Y,c g (s g ) + Y,c c (S c ) 

geg ceC 


(1) 


2 In the rest of the paper, complex numbers are denoted with a bar above (e.g., V) and 
complex conjugates with a bar below (e.g.,\f). 
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Fig. 1: Notation used in this paper for the OPF formulation. 


The first term of the network objective (C g ) in (1) is typically a non-decreasing 
convex function accounting for the minimization of the generation costs or the 
network real power losses. The second term ( C c ) is included in the objective 
when the cost of non-supplied load is taken into account. 

The following set of constraints is considered 3 : 
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3 Note that the proposed formulation can be extended without loss of generality to the 
case of multi-phase unbalanced grids by adopting the so-called compound network admittance 
matrix, (i.e., the 3-phase representation of the grid model which takes into account the various 
couplings between the network phases) instead of the single-phase equivalents. In this cases, 
each of the constraints in (2)-(9) needs to be formulated separately for each network phase. 
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where, S denotes the complex power 4 , Vb is the direct sequence phase-to-ground 
voltage of node b, I(+ (If -) is the current flow in the receiving (sending) end of 
line t, and 'Hg.'Hc are the capability curve of the generator g and the limits of 
the load c respectively 5 . If a generator (load) is non-controllable then the set 
T~L g (' H c ) is limited to a single point. 

The first constraint (2) corresponds to the power balance constraint at each 
network bus, whereas (3) is an alternative way to define the AC power flow 
equations. Constraints (6) and (7) are so-called node voltage and lines ampac¬ 
ity contraints, i.e., limits on node voltages and line power/current flows. The 
last constraints (9) represent the capability limits that each of the controllable 
devices should respect. 

The equality constraints (3) render the OPF problem non-convex and, there¬ 
fore, difficult to solve efficiently. The majority of the proposed algorithms in the 
literature rely on several approximations and/or convex relaxations and seek a 
solution to a modified OPF problem. In what follows, we describe and discuss 
the most common approximations. 

2.3. Approximations of the OPF Problem 

In general, the approximations used in the formulation of an OPF problem 
can be categorized in two large groups: approximations of the physical net¬ 
work models and methods that relax the space of the solutions and/or control 
variables. 

In the first case, we can find OPF formulations that rely mainly on lineariza¬ 
tions of the AC power flow equations. Such attempts typically (i) consider the 
DC power flow, (ii) use the decoupled AC power flow or (iii) neglect the net¬ 
work losses and/or the transverse parameters of the lines. Specifically, the 
concepts of the DC and the decoupled OPF have been extensively used in the 


4 We use the convention that positive values represent power injection and negative power 
consumption. 

5 Note that different types of controllable generators or loads can be accounted for via their 
corresponding capability curves/limits. 
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literature (e.g., [32, 33, 34, 35]), as they approximate the OPF problem with 
linear programming problems and, therefore, enable its fast resolution. Such 
techniques have been extensively used for the OPF solution in the case of trans¬ 
mission networks where the DC approximation might be reasonable since the 
resistance to inductance ratio of the transmission lines is negligible and the re¬ 
active power flow is supplied locally. However, the DC OPF is not applicable 
in distribution networks due to the large resistance over reactance ratio of the 
lines. Furthermore, the DC-based OPF algorithms have several shortcomings 
like, for instance, the inability to optimize the reactive power dispatch and the 
fact that they always provide a solution even when it is physically infeasible. 
In the same direction, the authors in [24] use the so-called Dist-Flow equations 
([36]) to linearize the power flows and propose an ADMM-based OPF algorithm 
that neglects the real and reactive losses. Finally, several contributions rely on 
simplified network line-models that neglect the transverse parameters, resulting 
in inaccuracies of the physical system model (e.g., [37, 38, 39]). 

In the second case, we can find OPF formulations where, typically, the con¬ 
straints are relaxed in order to convexity the problem. In particular, a large num¬ 
ber of contributions recently proposed a SDP formulation of the OPF problem, 
where the rank-one constraint of a matrix is relaxed and the algorithm is claimed 
to yield zero-duality gap for radial distribution networks (e.g., [40, 20, 21]). An¬ 
other relaxation is proposed in [39] where the OPF problem is cast as a second 
order cone programming. A similar technique is used in [41], where the equality 
constraints of the branch flows are relaxed. 

In both the aforementioned categories of approximations, the modified OPF 
formulations guarantee convergence of the proposed algorithms. The reached 
solutions, however, even though mathematically sound, are not always mean¬ 
ingful for the grid operation. The DC and the decoupled OPF work sufficiently 
well for transmission systems, nevertheless they can introduce large errors when 
used for solving the OPF in the case of distribution systems (e.g., [42]). As far 
as the semidefinite relaxation is concerned, its limitations have been recently 
investigated. The authors in [28] show through practical examples, that in the 



case of negative locational marginal prices or strict line-flow constraints it can 
lead to solutions that are not valid, namely for which the duality gap is not 
zero. Furthermore, in [43] the authors show the existence of multiple local op¬ 
tima of the OPF problem due to the feasible region being disconnected and due 
to the nonlinearities of the constraints; they show that the SDP formulation 
of the OPF problem fails to find the global optimum in cases where there are 
multiple local optima. In the same direction, a recent review ([38]) summarizes 
the semidefinite relaxations applied to the OPF problem and discusses their 
limitations. 

Recently, another formulation of the OPF problem has been proposed ([44, 
45, 29, 30, 46]). This formulation also belongs to the category of the semidefinite 
relaxations and uses the so-called branch-flow model (BFM) for describing the 
network. The BFM essentially describes the network flows by using as variables 
the currents and the powers of the various network branches, instead of the 
nodal injections. In [29, 30] Farivar and Low propose an OPF formulation that 
relies on the BFM representation of the network and they present a two-step 
relaxation procedure that turns the problem into a second-order cone program 
(SOCP). The authors prove that under specific assumptions both relaxation 
steps are exact for the case of radial networks, hence a globally optimal OPF 
solution can be retrieved by solving the relaxed convex problem. 

In what follows, we first briefly recall the formulation of the OPF problem 
in [29, 30] and then we investigate the applicability of the branch flow model 
to the OPF formulation. We show, on one hand, that this model misinterprets 
the physical network representation by imposing the ampacity constraint on a 
fictitious line-current that neglects the contribution of the shunt components of 
the line. We also show that, on the other hand, the proof of the exactness of 
the proposed relaxation requires several unrealistic assumptions. In particular, 
the OPF formulation in [29, 30] assumes full controllability of both loads and 
generators connected in the network buses. As a matter of fact, this is a strong 
assumption, as in a real setting the DNO has very few specific control points 
available in the network with controllable resources’ capability curves that are 
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typically complex. In addition to this, the controllable loads are required to 
have infinite upper bounds in order to prove the exactness of the proposed 
SOCP relaxation in [29, 30]. In a realistic scenario, such an assumption implies 
that in cases where excessive production of the generators causes violations of 
the voltage or line-flows limits, local demand is invoked to compensate for the 
increased generation even beyond the possible power consumption from installed 
loads. 

3. On the Limits of the Branch-Flow Convexification for the Solution 

of the OPF Problem 

3.1. The BFM-based Formulation of the OPF Problem 

We assume the same objective function as in Eq. 1 and again consider that 
the network lines are represented using a 7r-model. Contrary to the formulation 
in (2)-(9), we reformulate the constraints of the OPF problem by using the 
branch power and current flows as variables, similarly to [29]. To this end, we 
denote by S/ and If the power and the current that flow across the longitudinal 
elements of a network line I from the receiving toward the sending end, for 
which it holds that 


h = Yt(V pm -V p{e - ) ), W eC 

( 10 ) 

St = V/3(t+)lti Vf € C 

( 11 ) 


The power and current flows along the shunt elements of the lines are taken 
into account in the bus power balance constraints as nodal injections. In this 
direction, we denote by Yj, 0 the sum of all the shunt elements of the lines that are 
adjacent to bus b. In particular, the notation used for the BFM convexification 
of the OPF problem is shown in Fig. 2. 

Using this nomenclature, the constraints of the OPF problem are reformu¬ 
lated as follows: 
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Fig. 2: Notation used in this paper for the BFM convexification of the OPF problem. 
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Note that in this formulation of the OPF problem, the capability curves of 
the controllable loads and generators, i.e., constraints (17,18) on the nodal power 
S are limited to rectangular regions. This is essential for the conic relaxation 
proposed in [29, 30]. 

Starting from this formulation, the equality constraints in (13) are relaxed 
to inequalities and the aforementioned problem is casted as a second-order cone 
program. They also prove that for radial networks a global solution of the 
original OPF problem can be recovered from the solution of the relaxed problem 
if there are no upper bounds on the loads. In other words, the OPF problem is 
solved (12)-(18) by setting P Cmao; =oo and Q Crnaa; = oo in constraint (18). 

We show, in what follows, that this formulation is not equivalent to (1-9). 
In particular, constraint (16) (constraint (9) in [29]) is only an approximation of 
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the ampacity constraints and, moreover, the assumptions on the controllability 
and bounds of the energy resources in the network are unrealistic. 

3.2. Misinterpretation of the Physical Network Model in the BFM-based OFF 
Formulation 

The branch-flow model has been often used in load-flow studies (e.g., [47, 48]) 
and constitutes an accurate representation of the network model. The first 
problem with the formulation in (12)-(18) is that it misinterprets the physical 
network model when constraining the line flows in the network. Even though 
the power-flow equations in (12)-(14) are exact when the shunt capacitances 
are considered as nodal injections, the constraint (16) is imposed on a fictitious 
current flow across the longitudinal component of the lines, thus does not ac¬ 
count for the current flow toward the shunt elements. Therefore, the optimum 
of problem (12)-(18) can be such that the line ampacity constraint is violated. 

To better clarify why this occurs, we use a single-branch toy network, as 
shown in Fig. 3. The line parameters, as well as the base values of the system 
are given in Table 1. A purely resistive load is connected to bus 2 that we vary 
linearly in the range of [100 — 10000] Ohms in order to numerically quantify the 
mismatch between those quantities. We measure the current flows at the two 
ends of the line, as well as the flow along the longitudinal impedance of the 
line. Fig. 4 shows the measured quantities as a function of the load. It can be 
observed that the current flowing across the longitudinal impedance of the line 
under-estimates the actual current flow in the receiving end of the line. 

Bus 1 j Bus 2 

LF -1—* I-•lf 

I V | PI1 5 | Controllable_Load 


Fig. 3: The test network used for the numerical comparison of the current flows at the send¬ 
ing/receiving end of the lines and the current flow along the longitudinal line impedance. 
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Table 1: Parameters of the test network in Fig.3 


Parameter 

Value 

Network rated voltage, V (kV) 

15 

Line parameters, i?(Ohms), L(H), C(uF) 

(1,0.003,0.54) 



Fig. 4: Current flows at the sending/receiving end of the line and along the longitudinal line 
impedance (log-log scale). 



Controllable_device_1 Controllable_device_2 
Fig. 5: Network used in the study of the BFM-based OPF formulation. 


As a consequence, in this formulation of the OPF problem setting the limit 
on the longitudinal current flow below the line ampacity does not guarantee 
that the actual line current will respect this limit. In order to illustrate such 
a scenario, we consider yet another simple test network shown in Fig. 5. All 
the network lines are built by using the same values of resistance, reactance 
and capacitance per km, but by assuming different values of their length 6 . We 


6 Typical values of medium-voltage underground cables are considered for the resistance, 
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Table 2: Parameters of the test network in Fig.5 used for the investigation of the line ampacity 
limit violation_ 


Parameter 

Value 

Network rated voltage and base power, V(kV),S(MVA) 

24.9,5 

Line parameters, i?(Ohms/km), L(mH/km), C(uF/km) 

(0.193,0.38,0.24) 

Lines length, (km) 

(2.5,3,3.5) 

[P Smin ,Pg m J (MW) 

[0,2] 

P Cmin (MW) (bus2, bus3) 

(0.05,0.06) 

Q Cm in (Mvar) (bus2, bus3) 

(0.03,0.027) 

[^ram? ^raatc] (P-u) 

[0.9,1.1] 

Imax (A) 

80 


assume a first test case where the controllable device connected to bus 4 is 
a generator, whereas controllable loads are connected to buses 2 and 3. The 
network characteristics, the base values, the capability limits of the controllable 
resources 7 , and the voltage and ampacity bounds are provided in Table 7. We 
assume that the controllable generation operates at a unity power factor. The 
problem in (12)-(18) is formulated and solved in Matlab. The objective function 
accounts for loss minimization, as well as utility maximization of the controllable 
generation units: 


min - y^Re^Sg) + Re(Y e )\I e \ 2 (19) 

s 9 ’ s *’\v t \,\u\ 

In order to investigate the order of magnitude of the violation of the ampacity 
constraint, we solve the OPF problem for various line lengths and network 
voltage-rated values. In particular, we assume that the line lengths are uniformly 
multiplied by a factor in the range [1.25 — 7.5] (while keeping the network voltage 


reactance and shunt capacitances of the lines taken from [49], 

7 The upper bounds of the active and reactive power of the loads are considered to be 

infinite, as required in [29, 30]. 
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Table 3: Parameters of the test network in Fig.5 used for the investigation of the network 
operating point on the line ampacity limit violation 


Parameter 

Value 

I4„.4J(MW) (bus 2) 

[0,0.01] 

[Pg m ^Pg m J(MW) (bus 3) 

[0,0.012] 

(P Cmi „,Q Cmi J(MW,Mvar) (bus 4) 

0.3,0.15 


rated value to its nominal value) and the network voltage rated value varies in 
the range [15 — 45]AV (while keeping the line lengths to their nominal values). 
Once the optimal solution is computed in each case, we calculate the actual 
current flows in the sending/receiving end of the lines and we compute the 
maximum constraint violation. The results are shown in Fig. 6. As the line 
length increases, the current flowing toward the shunt capacitors increases, thus 
neglecting its contribution to the line flow leads to significant violations of the 
ampacity limit. At 7.5 times the initial line length, the violation reaches a 
value of 18.4%. The effect of the network voltage-rated value is similar, with a 
maximum constraint violation of 25% when the voltage value is 45kV. 

In addition to the effect of the line lengths and the network voltage-rated 
value, we study the effect of the network operating point on the ampacity viola¬ 
tion. To this end, we consider a second test case where the controllable device 
connected to bus 4 is a load and generators are connected to buses 2 and 3. 
The capability limits of the controllable resources are provided in Table 3. For 
this setting, Fig. 7 shows the solution of the BFM-based OPF problem, namely 
current flows at the receiving/sending end of the network lines, as well as across 
the longitudinal impedance. We can observe that the maximum violation of the 
ampacity constraint is in the order of 39.6%. 

In order to avoid current flows that exceed the lines’ ampacity limits, i.e., in 
order to use the BFM in an accurate way, the aforementioned formulation should 
either consider the actual current flows in the receiving/sending ends of the lines 
as optimization variables, or should add the contribution of the current flows 
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Line length factor Network voltage rated value (kV) 

Fig. 6: Maximum ampacity constraint violation as a function of the line lengths and the 
network voltage rated value. 



Line index 


Fig. 7: OPF solution for the current flows at the sending/receiving end of the network lines 
and across the longitudinal line impedance under heavy consumption and light generation 
conditions. 
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toward the shunt elements of the lines to the longitudinal current flow in the 
inequality constraint (16). By adopting either of the two approaches, however, 
(12)-(18) can no longer be solved efficiently as proposed in [29, 30]. Therefore, 
the generic OPF problem cannot be convexified by using the approach in [29, 30]. 

3.3. On the Assumptions Required for the Exactness of the SOCP Relaxation 

In addition to the aforementioned fundamental problem, which is related 
to the physical network model, the proof of exactness of the proposed SOCP 
relaxation in [29, 30] requires specific assumptions, related to the controllability 
of the demand in the network. Several of these assumptions might not be 
realistic. The goal of this section is to discuss these assumptions and their 
consequences in a realistic setting. In order to do so, we consider realistic 
case-studies and we show that the solution of the OPF problem can result in 
unrealistic values for the control variables. 

To begin with, the OPF formulation in [29, 30] assumes controllability of 
both loads and generators in the network buses and, in particular, assumes 
rectangular bounds on the powers of loads/generators. This is quite a strong 
assumption, as usually the DNO has very few specific control points available 
in the network with capability curves that are typically more complex and that 
account, among others, for capabilities of power electronics and limitations of 
machinery. An even more serious limitation is that the model in [29, 30] consid¬ 
ers no upper bounds on the controllable loads in order to prove the exactness 
of the proposed relaxation. This implies that in cases where excessive produc¬ 
tion of the generators causes violations of the voltage or line-flows limits, local 
demand is invoked to compensate for the increased generation. In order to il¬ 
lustrate such a setting and to show that the result of the OPF problem can 
result in unrealistic values for demand, we consider the same network in Fig. 5 
and we assume that there is high penetration of distributed generation and a 
low demand. The values of loads and generation, as well as the corresponding 
limits are shown in Table 4. Solving the optimization problem and considering 
infinite upper bounds on the demand results in load values that are significantly 
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Table 4: Parameters of the test network in Fig.5 used for the investigation of the unbounded- 
ness of the consumption 


Parameter 

Value 

[P gmtn ,P gm J (MW) 

[0,1.2] 

(Pc min ,Qc min ) (MW,Mvar) (buses 2,3) 

(0.0125,0.0026) 


increased, compared to the minimum values shown in Table 4. The resulting 
optimal power points are shown in Fig. 8. We show in black the initial values for 
active and reactive power of loads and generation (corresponding to the values 
of Table 4), and in gray the results of the OPF solution (when not accounting 
for upper bounds on loads). It is worth observing that the optimal active power 
consumption of bus 3 is increased 23.6 times and the reactive power consump¬ 
tion at buses 2 and 3 is increased 85.3 and 92 times, respectively. In a realistic 
setting, even if part of the demand in the network is controllable, the amount 
of available demand-response is limited and such an increase in the consump¬ 
tion is most likely not possible. Therefore, in such a case, the congestion and 
voltage problems should be solved by properly controlling the generator within 
its capability limits. In addition to this, typically, the active and reactive power 
consumption should be linked via the corresponding power factor. We observe, 
however, that the OPF solution in this scenario results in very large values for 
the reactive power consumption and, in particular, the power factor of bus 2 is 
0.03 after the OPF solution, whereas initially its value is 0.98. In an attempt to 
relax this assumption, it is shown in [46] that the infinite upper bound on the 
loads, when not applicable, can be replaced by equivalent conditions. However, 
not only are these conditions unrealistic, they are also not applicable in our 
context as they require no upper bound on the voltage magnitudes. This is in 
contradiction with the actual problem we target, i.e., voltage rise due to high 
penetration of renewable energy resources. 


18 





Bus index 


x 10 


' ' ^■initial minimum (maximum) value 
loads (generation) 

6 - l lf M Farivar-Low OPF solution with 
infinite upper bounds on loads 


I 4 


■S3 


1 2 3 

Bus index 


Fig. 8: Optimal solution of the OPF formulation for the active and reactive power set-points 
when upper bounds on loads are infinite. 


3.f. On the Extension of the SOCP Relaxation to Networks with Lines Modeled 
as 7r- equivalents 

In this paragraph, we discuss two different approaches that can be used to 
extend the initial formulation in (12)-(18) in order to properly account for the 
shunt elements of the lines and the line ampacity constraints. The goal of this 
paragraph is to show through concrete examples that extending the approach 
in [29, 30] to a system model that is correctly represented results in a convex 
problem, however it cannot guarantee the exactness of the SOCP relaxation 
and, thus, the retrieval of a feasible OPF solution. 

The first straightforward way to properly account for the line ampacity con¬ 
straints in (16) is to include in these inequality constraints the contribution of 
the current flowing towards the shunt elements of the line. In order to do so, 
we keep the same branch flow variables of the initial formulation in (12)-(18) 
and we define a new set of line constraints for the case of 7r-model lines. In this 
case, the total line flowing, for instance, in the receiving end of the line is the 
sum of the longitudinal current plus the contribution of the shunt, resulting in 
the following constraint: 

|h + Y tt V w+) | 2 < 7? mox , W € C (20) 
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Similarly for the sending end of the lines. Now, after expanding the square, 
the line constraints (16) are reformulated as follows 8 : 

\h \ 2 + \Y e +\ 2 \Vp (e+) \ 2 + 2Re(Y e +S e ) < W G C (21) 

\h\ 2 + \Y t - | 2 |C W -)| 2 + 2Re(Y e -(Yf 1 \Ie\ 2 - S e )) < lf max , W €E £ (22) 

It is worth noting that these new constraints, that account also for the line 
flows towards the shunt elements, are still convex in the branch-flow variables, 
i.e., (|Vs(^+)| 2 , \Ie\ 2 , Re(Se), Im(Se)). Therefore, these constraints can be added 
to the initial formulation in [29, 30] without losing convexity. The main question 
now is whether the SOCP relaxation continues to be exact with the reformula¬ 
tion of the lines ampacity constraints as above, i.e., whether the optimal solution 
of the relaxed SOCP problem is guaranteed to be a physically feasible one. As 
we show below, there are cases for which the attained solutions can be physi¬ 
cally infeasible. In other words, even with a correct system model, the proposed 
SOCP relaxation in [29, 30] cannot guarantee a solution that is meaningful for 
the grid operation. 

Let us consider the following simple example. We use once again the simple 
test network shown in Fig. 5. We assume a test case where the controllable de¬ 
vice connected to bus 4 is a load, whereas controllable generators are connected 
to buses 2 and 3. The network characteristics, the base values, the capability 
limits of the controllable resources, and the voltage and ampacity bounds are 
provided in Table 5. We assume that the controllable generators operate at a 
unity power factor. Note that the upper bounds for the loads are considered 
infinite as required in the original formulation in [29, 30]. 

We solve the problem in Matlab, using the interior-point algorithm provided 
by the fmincon solver, and the resulting values for the SOCP inequalities for 


* Note that with the inclusion of the shunt elements of the lines two inequality constraints 
are required per line in order to properly account for the ampacity limits. The reason is that 
in this case, the currents at the two ends of the line are no longer equal and both need to be 
constrained below the line ampacity limit. 
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each network line are shown in Table 6. In this case, after the solution of the 
OPF problem not all the inequalities in (13) are satisfied with equality, namely 
the SOCP relaxation is inexact, therefore, the obtained solution has no physical 
meaning and a physically feasible solution cannot be recovered. 

A second way to extend the approach in [29] to networks with lines repre¬ 
sented as 7r-equivalents is to reconstruct the BFM in order to include the shunt 
elements of the lines. To this end, we consider an undirected graph G = (Af,£). 
A f and 8 represent the set of nodes and lines respectively. The power flow equa¬ 
tions that define the branch-flow model with the inclusion of shunt elements are 
in this case: 


E 

4-E^ 

-E 5 * + 

-£*<- = 

0, V6 G B 

(23) 

g£b 



P(l~)=b 



s i+ 

= Vp (e+)le+ 

, S e - = 


\/£gC 

(24) 


= 

1 

1 

+ Y t +Vp {l +), 

\/£eC 

(25) 

h- 

= YtiVw-) 

(Ns 

+ 


\/£gC 

(26) 


Note that, contrary to the formulation in(10)-(ll), the power flows and 
current flows variables are now defined for both ends of the lines. 

Similarly to the procedure described in Section III.A in [29], we substitute 


Table 5: Parameters of the test network in Fig.5 


Parameter 

Value 

Network rated voltage and base power, P(kV),£(MVA) 

24.9,5 

Line parameters, iZ(Ohms/km), L(mH/km), C(uF/km) 

(0.193,0.38,0.24) 

P gmax (bus2, bus3) (MW) 

(1,1.2) 

P Cmin (MW) (bus4) 

0.1 

Q Cmi „ (Mvar) (bus4) 

0.05 

[^rain? ^raatc] (p*u) 

[0.9,1.1] 

1-max (A) 

80 
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Table 6: SOCP inequalities in (13) 


Line 

Value 

1-2 

1.7 E- 16 

2-3 

-0.0461 

3-4 

1.3 E- 16 


(25) and (26) in (24) resulting in the following set of equations per line: 

(Ye + %)^ 2 (/+ ) = S e+ + YtVpwVpy-) (27) 

(Ye + Y t - )Vfo~) = S t - + Y e V m) V^-) (28) 

Then taking the magnitude of the resulting equations squared, we can elim¬ 
inate the angles from (27)-(28). Doing so, the OPF problem in (12)-(18) is 
reformulated as: 


_ min VW + E^) 

S a ,S c ,S < +,|I r |.| ,|Vg (f+ j] , ge g cgC 

S e - ,| If- | 2 iV/3(<-)| 2 


(29) 

subject to: ^ S g — ^ S c — ^ Sg+ — ^ Sg- = 0, V6 € B 

gGb cGb /3(£+)=6 0 (e~)=b 


(30) 

\Sg+ 2 - \Vp(e+)\ 2 \h+\ 2 < 0, \Sg-\ 2 - \ V p ( t -)\ 2 \ T t -\ 2 < 0, 

V£e £ 

(31) 

K +| 2 | V ^ +) | 2 - \V K g-)\ 2 = 2 Re(a e+ Yj 1 Sg + ) - \Y^\\Ig + \ 2 , 

WieC 

(32) 

| a /-| 2 |^-)| 2 - \Vp { g +) \ 2 = 2Re(a e -Yj 1 Sg-) - | YT 1 ||/,-| 2 

, C 

(33) 


where ag+ := 1 + Y^ _1 Y^t-, ag- := 1 + Yf^Y^. 

At this point it is important to note that, in order to recover a solution of the 
original ACOPF problem, we need to recover the line angle from the solution of 
the above problem in a way similar to [29]. In the original paper, it is shown that 
for radial distribution networks the angle relaxation step is always exact. On the 
contrary, in the formulation above, both angles f3g+ = Z(ag+\Vp/g+\\ 2 — Yj 1 Sg+) 
and fig- = Z(gg-\Vp(£-)l 2 — Yj 1 Sg~) are defined. In order for a solution to be 
physically meaningful, namely the angle relaxation step to be exact, these line 
angles should satisfy (3g+ + fit- = 0. However, there is no guarantee that this 
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Table 7: Parameters of the test network in Fig.5 


Parameter 

Value 

Network rated voltage and base power, P(kV),S(MVA) 

24.9,5 

Line parameters, i?(Ohms/km), L(mH/km), C(uF/km) 

(0.193,0.38,0.24) 

[P gm ^P 9m J (MW) 

[0,2] 

P C „„(MW) (bus2, bus3) 

(0.05,0.06) 

Q Cmin (Mvar) (bus2, bus3) 

(0.03,0.027) 

[^rain? ^raatc] (P-U) 

[0.9,1.1] 

Imax (A) 

80 


will occur in the obtained solution. In fact, as we show in the example that 
follows the angle relaxation is not exact when using this formulation even in the 
case of radial networks. 

To support the above claim, we consider the same simple test network shown 
in Fig. 5. For the sake of simplicity, we assume a test case where the control¬ 
lable device connected to bus 4 is a generator, whereas controllable loads are 
connected to buses 2 and 3. The network characteristics, the base values, the 
capability limits of the controllable resources, and the voltage and ampacity 
bounds are provided in Table 7. We assume that the controllable generation 
operates at a unity power factor. Note that the upper bounds for the loads are 
considered infinite as in the original paper ([29]). 

The objective function has two terms, namely loss minimization and utility 
maximization of the controllable generation unit: 

_min ^f?e(F r 1 )(|/ f+ | 2 + \V e+ \ 2 V/+ - Im(S e+ ))+ (34) 

S g ,S c ,S e+ ,\I e+ \ 2 ,\V e+ \ 2 ,f-t 

S e -,\I e -\*,\V e -f 

ReiYf'XlIe-l 2 + \V e -\% 2 - - Im(S ( -)) -Y,Re(S g ) (35) 

g£Q 

It is worth mentioning that in view of the new formulation, the current used 
for the real losses computation is no longer represented by the variable 7? | 2 , but 
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needs to be computed as the difference between the current flowing from one 
end of the line to the other and the current flowing towards the shunt elements. 
This has two implications. First, the current across the series admittance can be 
computed twice using | Ig+ 1 2 or \Ig- | 2 . Neglecting one of the two currents in the 
objective function results in a non-exact SOCP relaxation. Second, computing 
the longitudinal component of the current results in the objective function not 
being independent of the power flow variables, which is one of the assumptions 
used in [29, 30] to prove exactness of the proposed relaxation. 

In this case, after the solution of the OPF problem the inequalities in (31) 
are satisfied with equality, namely the SOCP relaxation is exact. The values for 
the two inequalities for each network line are shown in Fig. 9. For the same test 
case, we obtain the line angles /3^+ and /3^-. The results are shown in Fig. 10. 
One can clearly observe that the obtained line angles do not satisfy /3^+ +0t- = 0 
and, therefore, the obtained solution has, again, no physical meaning. 

The aforementioned examples indicate that the proposed SOCP relaxation 
in [29, 30] cannot be trivially extended to lines represented as 7r-model equiva¬ 
lents even with convex constraints applied to the line exact 7r-model. 



Line index 

Fig. 9: The relaxed inequalities values for the problem (29)-(33). 
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Line index 


Fig. 10: The line angles obtained after the OPF solution of problem (29)-(33). 

3.5. Discussion 

In the previous sections, we have investigated the BFM-based OPF formu¬ 
lation that allows a SOCP convexification of the OPF problem, therefore an 
efficient resolution of the problem with guaranteed convergence. This particu¬ 
lar convexification proposed originally in [29, 30] is valid in specific cases. In 
particular, when the DNO knows a-priori that the network lines are operated 
far away from their ampacity limit, therefore constraints (16) are irrelevant. 
Finally, when the network lines are electrically short, namely when their shunt 
elements, combined with the line operating voltage, drain a negligible amount 
of capacitive reactive power. Practical examples are composed by short over¬ 
head lines operating at nominal voltages less than 20kV. However, as we have 
shown in the previous sections there are several case-studies involving realistic 
network topologies and operating points for which the original formulation pre¬ 
sented in [29, 30] results in violations of the line ampacity limits. Additionally, 
the assumptions required for the exactness of the proposed relaxation can re¬ 
sult in solutions that are unrealistic for the grid operation, when for instance 
increased connection of distributed generation units are connected to the grid. 
Finally, we have extended the original problem formulation to include the shunt 
capacitances in an effort to properly account for their contribution in the line 
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ampacity constraints of the OPF problem. As we have shown, such an ex¬ 
tension preserves convexity, however cannot guarantee exactness of the SOCP 
relaxation. In this work, we are interested in the generic formulation of the 
OPF problem without any restriction on the grid topology and its operating 
point. As a consequence, there is a need to design algorithms that target the 
original non-approximated OPF problem that remains inherently non-convex. 
Recent trends are in favor of using ADMM for the solution of the OPF problem. 
Even though ADMM requires the underlying problem to be convex in order to 
guarantee convergence, it has been applied also to the case of non-convex AC 
OPF problems with promising convergence performance (e.g., [23, 26]). In what 
follows we first present the ADMM solution of the problem in (l)-(9) and then 
we highlight specific scenarios for which ADMM fails to converge when applied 
to the non-approximated OPF problem. 

4. On the Application of ADMM for the Solution of the OPF Problem 

f.l. ADMM-based Solution of the OPF Problem 

The ADMM-based solution of the OPF problem requires that the control 
variables are split into two separate groups and that the objective function is 
separable across this splitting [50]. To this end, we introduce additional slack 
variables, z, for the devices’ and loads’ power injections and for the line power 
flows and we reformulate the OPF problem as follows 9 : 


9 In what follows we assume that demand is non-controllable. Also, as in [22] the con¬ 
straints (3),(9) are considered internal constraints of the lines and devices respectively and 
if , 17 are internal variables of the lines. 
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min -£tf fl (ife(S fl )) + £jy(|H|)+ (36) 

Sg,Zg^S c ,Zc ,S £ +,Z £ + ,S t - , 

z £ - ,E £ + ,E e - ,I e + ,I t - Vb 

52 < MI51> 151)+£ $(52 *9 - 52 * c +52 *<■++£ **-) 

^ b g&b ceb /3(£+)=b fi{t.~)—b 

subject to: S g = z g ,\/ g £ Q, and S c = z c ,\/ c £ C (37) 

= Zg +, and 5^- = Zg_- ,\/i£C (38) 

Ee+ = Vp(e+)f an< i ^e- = Vp{t-)i XI l £ C (39) 

where 0 is the characteristic function of the set {ieC:i = 0}, Jy is a penalty 
function with value 0 if V m i n < | V), < 17 maa . and oo otherwise and Jj is a penalty 
function with value 0 if max(|/^“|, |7jT|) < Ie max and oo otherwise. 

The augmented Lagrangian for this problem is as follows: 

S c , Sg+ , *S^— , Eg+ , , Ig+ , Ip- . Z g , Z c , 2^>+ , 2^— , Vb, p., b', A) 

= - £ U s{MS g )) + £ Jy(|V0|) + £ J/(|7/|, I5“l) 

9 b l 

+ £ <£(£ % - £ *= + £ £«/-) 

b geb cSb /3(€+)=b P(l~)=b 

+ ^{£ 1-^5+ - ^9(*+) + /£| 2 + £ |£*- - Vg(f-) + Vr| 2 
£ £ 

+ £ l^g — Zg + ^g | 2 + £ \&c ~ Z c + A c | 2 

g c 

+ |«S^+ — zp+ + A^+1 2 + |5^+ — Zf- + Af- |"} (40) 

e e 

where Jl. v 1 A are the lagrange multipliers associated with the equality constraints 
(37)-(39). 

The ADMM algorithm at the k—th iteration consists of the following steps: 

1. First, all the devices, loads and lines update in parallel the primary vari¬ 
ables, and their internal variables, i.e., ( S g , S c , Sp+ , Sp- , Ep+ , Ep- , h+M 
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with the secondary variables, and the dual variables fixed 10 


For each network line b. 

(qk+1 nk+1 fpk+1 fpk -\-1 fk +1 jk-\- 1 \ _ 

\°£+ ? 7 ■ Cj £+ i - C/ £- 5 1 f>+ 5 J g- ) ~ 

argrnin _ J/(|7+|, |/7|)+ 

(- — ,£/£_)_ — j I £-\-1 1 g — 

^ {\Ei+ - Vp(e+) + fiel 2 + \Ee- ~ Vp(i-) + I 2 


+ |5 /+ - 4- + X(+ 1 2 + | S t - - z\- + Af_ | 2 ) (41) 

subject to: = E(+I e + and St- = (42) 

h+ = Y e (E e+ -E t -) + Y ( + E e+ (43) 

h- = %(E e - -E t+ ) + Y t -E t - (44) 

For each device g: (45) 


Sg +1 = argmin - U g (Re(S g )) + ^(|5 S - z k g + A*| 2 ) 

s g z 

subject to: S g G T~L g 

For each load c: £k +1 = S c (46) 

2. Then, by using the updated primary variables, the secondary variables are 
updated, i .e.,(z, V&), on a bus level. We denote by Zi, the vector of complex 
powers of all the devices, loads and lines that are connected to bus b, i.e., 

Zb — {z g : g £bi Z c:C £b> Zg+ : f3(g+')—b, Zg- : f}(e-')—b)'- 

Zb +1 = argmin(0(E z g - E z c + E z e + E z t -) (47) 

g&b cGb py+)—b py-j-b 

+ |{£ l^9 +1 - ^ + A^l 2 + E \Sc +1 ~ ^ + A^l 2 

g£b cG6 

+ E — zt+ + ^e+ 1 2 + E - 1 2 }) 

/3p+)=b /3 (t~)=b 


10 Note that demand is not controllable, hence the loads do not require the solution of an 
optimization problem to update their power consumption. 
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( 48 ) 


V b k+1 = argmin( J(V b ) + |^ fc + +1 - V„ + ^ 

Vb P(t + )=b 

+ £ I E k + l - V b + ^ fe | 2 }) 

0 (t~)=b 

3. Finally, dual variables, i.e., /2, P, A are updated: 

£/ +1 = jS? + (Efr 1 - V k +l } ) (49) 

P* +1 =^ + (^ fc _ +1 -V;+_ 1 ) ) (50) 

A g fc+1 = Ag + (S k+1 z k+1 ) (51) 

Ac +1 = Ag + (S k+1 z k+1 ) (52) 

A^ 1 = A^+ + (^ fc + +1 - 4+ 1 ) (53) 

A^ 1 = X k - + (5 k + x - z k + l ) (54) 

The stopping criterion for this algorithm is that the primal and dual residuals 
(defined as in [50]) are less than a small predefined tolerance or that a maximum 
number of iterations has been reached. 

In what follows, we show specific scenarios where the AD MM algorithm fails 
to converge to a solution. 

4-2. Investigation of the Convergence of the ADMM-based Solution of the OFF 
Problem 

We consider the same network in Fig. 5. Each network bus, apart from the 
slack bus, has a load and a generator connected to it. The demand in the net¬ 
work is assumed to be non-controllable, whereas the generators are assumed to 
be distributed solar panels with typical PV-type capability constraints. For this 
scenario, the capability limits and the values of loads and generation are given 
in Table 8. In addition to the loads and generation, we consider that a shunt 
capacitor is connected to bus 2. In order to model this shunt capacitor, we 
consider that it is part of the first network line. In particular, we consider that 
the shunt capacitance on the sending end of the 7r-model of the line connecting 
buses 1 and 2 is modified accordingly to account for the shunt capacitor. It is 
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worth noting that switched capacitor banks, if present in the network, can be 
taken into account in the problem formulation in a similar way provided they 
are not included in the OPF control variables. The control of these discrete ele¬ 
ments is possible provided that we assume they are continuous control variables 
which are rounded to the nearest integer upon solution of the OPF problem. 
In particular, switched capacitor banks, can be readily taken into account as 
control variables in the OPF formulation and in particular as nodal injections 
with controllable reactive power in a range defined by their capacity limits. 

We implement and solve the ADMM algorithm in Matlab for two different 
cases that correspond to two different values of the size of the shunt capacitor 
(see Table 8). In Case I, even though the OPF problem solved is the non- 
approximated non-convex one, ADMM converges, within the predefined toler¬ 
ance, in 411 iterations. The left figure in Fig. 11 shows the objective function 
value as a function of the number of iterations of ADMM. The left figure in 
Fig. 12 shows the convergence of the buses’ voltage magnitudes and Fig. 13 
shows how the primal and dual residuals evolve with the iterations. On the 
contrary, in Case II, ADMM fails to converge to a solution and reaches the 
maximum number of iterations. This is shown in Fig. 11 (right), 12 (right) 
and 14 where the objective function, as well as the residuals and bus voltages 


Table 8: Parameters of the test network in Fig.5 used for the ADMM-based solution of the 
OPF problem 


Parameter 

Value 

Generators’ power, \Si gmax |, i = 2,3,4 (MVA) 

0.40,0.39,0.46 

Generators’ power factor, cos4>i g ,i = 2,3,4 

0.9 

Loads’ active power, P, c ,* = 2,3,4 (MW) 

2.76,2.16,2.46 

Loads’ reactive power, Qi 0 ,i = 2,3,4 (MW) 

1.38,1.08,1.23 

Shunt capacitor (bus 2), case I and II (uF) 

(239,859) 

Penalty term gain, u; 

1 

Tolerance and maximum number of iterations 

10~ 4 ,10 4 


30 






Fig. 11: Objective function value for case I and II (last 500 iterations). 
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Fig. 12: Voltage magnitude evolution for cases I and II (last 500 iterations). 


are plotted for the last five hundred iterations until the maximum number of 
iterations is reached; we can observe that they exhibit oscillations. 

In what follows we analyze why the ADMM algorithm converges in Case I 
but fails in Case II. To begin with, the first network line has the peculiarity that 
the voltage at its receiving end E(+ (i.e., the slack bus voltage) is fixed. 11 As a 
consequence, the first equality constraint in (42) becomes linear in the real and 
imaginary part of the voltage E ^~, whereas the second equality constraint in 
(42) becomes quadratic on the real and imaginary part of the voltage E(-. In 
fact, the coefficients of the quadratic terms in the latter constraint are Re(Y () 
and —Im(Ye) — Im(Y e ~) for the real and imaginary parts, respectively. Due 
to the physics of the network, Re(Yg ) and Im{Y ( ~) are positive for a network 
line and /m(Y)) is negative. Furthermore, typically, the longitudinal reactance 
Im{Yi) is much larger than the shunt capacitance Im(Y f -) of a network line. 


1 1 This holds for all the lines that are connected to the slack bus. 
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Fig. 13: Norm of the primal/dual residuals for case I (last 311 iterations). 
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Fig. 14: Norm of the primal/dual residuals for case II (last 500 iterations). 


Therefore, typically the coefficients of both quadratic terms are positive, and 
the line problem in (41) is convex for the lines that are connected to the slack 
bus. This is the case for the Case I. However, in Case II the size of the shunt 
capacitor, connected to bus 2, is such that Im(Y e ~) > —/m(Y^), thus the co¬ 
efficient of the aforementioned quadratic term in (42) is no longer positive and 
the corresponding line problem becomes non-convex. 

Apart from the aforementioned case of the shunt capacitor, we also discuss 
the case of OLTCs, which is another discrete control typically used to optimize 
the grid operation. The ADMM algorithm also fails to converge to a solution 
when on-load tap changers (OLTCs) are included in the OPF formulation as 
control variables 12 . To better understand why this occurs, let us consider a 


12 For the sake of brevity we do not include the simulation results for this specific scenario. 
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transformer with OLTC capabilities between buses 1 and 2 in the network and 
let us denote the ideal transformer admittance by Y t and the OLTC ratio by 
a £ R. Then based on the OLTC model in [51], the longitudinal admittance 
of the first network line equals aY t and the shunt elements of the receiving 
and sending ends of the same line are a(a — 1 )Y t and (1 — a)Y t respectively. 
If the OLTCs are not included in the set of control variables, then the ratio 
a has a fixed value and the inclusion of the OLTCs in the OPF formulation 
does not affect the solution. However, when the OLTCs are considered control 
variables 13 , their effect is similar to that of the shunt capacitors, in the sense 
that the line problem in (41) becomes once again non-convex for those lines 
that are connected to regulating transformers. The reason is that a is now 
an additional control variable, namely the OLTC ratio appears in the equality 
constraints (42) of the first network line problem, and both these constraints 
become quadratic in Eg- and a and non-convex. 

In this first part of the paper we have focused on investigating the limits of 
the branch flow convexification proposed by Farivar-Low in [29, 30] and of the 
ADMM-based solution of the OPF problem. In particular, we have discussed 
the misinterpretation of the physical model in the Farival-Low formulation of 
the OPF problem and the unrealistic assumptions therein. Finally, we have 
provided the ADMM-based decomposition of the OPF problem and we have 
shown, through specific examples, cases for which the ADMM-based solution of 
the non-relaxed OPF problem fails to converge. 
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Abstract 

In the first part of this two-part paper we show that the branch-flow convexifi- 
cation of the OPF problem is not exact and that the ADMM-based decompo¬ 
sition of the OPF fails to converge in specific scenarios. Therefore, there is a 
need to develop algorithms for the solution of the non-approximated OPF prob¬ 
lem that remains inherently non-convex. To overcome the limitations of recent 
approaches for the solution of the OPF problem, we propose in this paper, a 
specific algorithm for the solution of a non-approximated, non-convex AC OPF 
problem in radial distribution systems. It is based on the method of multipliers, 
as well as on a primal decomposition of the OPF problem. We provide a cen¬ 
tralized version, as well as a distributed asynchronous version of the algorithm. 
We show that the centralized OPF algorithm converges to a local minimum 
of the global OPF problem and that the distributed version of the algorithm 
converges to the same solution as the centralized one. Here, in this second part 
of the two-part paper, we provide the formulation of the proposed algorithm 
and we evaluate its performance by using both small-scale electrical networks, 
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as well as a modified IEEE 13-node test feeder. 

Keywords: OPF, ADMM, decomposition methods, method of multipliers, 
convex relaxation, active distribution networks, distributed algorithms, 
asynchronous algorithms. 


1. Introduction 

In Part I of this two-part paper we present the generic formulation of the 
non-convex OPF problem and we briefly review several OPF algorithms that are 
based on approximations and assumptions in order to guarantee convergence. 
Furthermore, we focus on the branch-flow convexification of the OPF problem 
that has been recently proposed by Farivar and Low in [1, 2] and is claimed to 
be exact for the case of radial distribution systems under specific assumptions, 
despite the absence of apparent approximations. We show that this claim, in 
fact, does not hold, as it leads to an incorrect system model and therefore, there 
is a need to develop algorithms for the solution of the non-approximated OPF 
problem that remains inherently non-convex. In detail, we show through prac¬ 
tical examples that in [1, 2], on one hand, there is a misinterpretation of the 
physical network model related to the ampacity constraint of the lines’ current 
flows and, on the other hand, the proof of the exactness of the proposed relax¬ 
ation requires unrealistic assumptions related to the unboundedness of specific 
control variables. Furthermore, we investigate the application of ADMM for the 
solution of the original non-approximated OPF problem. Even though ADMM 
requires the underlying problem to be convex in order to guarantee convergence, 
it was applied also to the case of non-convex AC OPF problems with promising 
convergence performance (e.g., [3, 4]). However, we show, through practical 
examples, cases for which the ADMM-based decomposition of the non-relaxed 
OPF problem fails to converge. 

To overcome the aforementioned limitations, here in this second part, we pro¬ 
pose an algorithm for the solution of the non-approximated non-convex AC OPF 
problem in radial networks. A large number of algorithms can be found in the 
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literature that tackle the non-approximated non-convex OPF problem ranging 
from non-linear and quadratic programming techniques, Newton-based meth¬ 
ods, interior point methods to heuristic approaches based on genetic algorithms, 
evolutionary programming, and particle-swarm optimization (e.g., [5, 6, 7, 8]). 
These methods, typically utilize powerful general purpose solvers or in-house 
developed software but they do not provide, in general, convergence guarantees. 
Our proposed solution belongs to the family of augmented Lagrangian methods 
for the solution of the OPF problem. More specifically, our proposed method 
uses an augmented Lagrangian approach, relies on the method of multipliers 
([9, 10, 11]) and provides convergence guarantees. In particular, we design a 
centralized OPF algorithm that is proven to converge to a local minimum of the 
original non-approximated OPF problem. 

With respect to the case of controlling multiple dispersed energy resources, 
it is of interest to also define a distributed solution method that is formally 
equivalent to the centralized formulation. Distributed solutions are of interest 
in several practical cases. Among others, when the problem size is large due to 
a very large number of small controllable resources like PV panels, when the 
communication requirements are such that the amount of exchanged information 
flow between agents needs to be limited, or when an asynchronous solution of 
the problem is more appealing. In fact, several distributed OPF algorithms 
are proposed in the literature. In [12, 13] the authors design a dual-ascent 
algorithm for optimal reactive power flow with power and voltage constraints. 
In [14, 15] dual decomposition is used as the basis for the distributed solution of 
the OPF problem. A further category of distributed OPF has been proposed to 
solve multi-area objectives ([16, 17]). However, their applicability to generically 
decomposable OPF has not been discussed. Finally, a significant number of 
contributions propose distributed formulations of the OPF problem, based on 
the alternating direction method of multipliers (ADMM) (e.g., [18, 14, 3, 19, 
20, 4]). 

In this direction, we present, here in this second part, a distributed version 
of the proposed algorithm that, unlike ADMM, is based on a primal decomposi- 
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tion [21] and does not require that the problem be convex. In this decentralized 
version of the algorithm, at each iteration, local agents, assigned to network 
buses and network lines, exchange messages with their neighbors using only lo¬ 
cal information. We prove that the distributed algorithm converges to the same 
solution as the centralized version. Finally, we present an asynchronous imple¬ 
mentation of the distributed algorithm where the messages of the neighboring 
agents need not be synchronized. 

The structure of this second part is the following. In Section II we describe 
the proposed algorithm for the OPF solution. We present both a centralized, 
as well as a decentralized asynchronous version of the proposed algorithm. In 
Section III we investigate the convergence of the proposed algorithm in the 
cases where the BFM convexification leads to an incorrect solution and ADMM 
fails to converge to a solution. In Section IV we evaluate the performance of 
the proposed algorithm using a modified IEEE 13-node test feeder. Finally, in 
Section V we provide the main observations and concluding remarks for this 
Part II. 

2. AC OPF in Radial Distribution Systems 

We first write the AC OPF problem presented in Part I in an equivalent 
form, and then we provide a centralized, as well as a distributed algorithm for 
its resolution. 

We make the following assumptions about the grid model: 

Al. We consider a direct sequence representation of the grid 1 ; 


1 Note that the proposed formulation can be extended without loss of generality to the 
case of multi-phase unbalanced grids by adopting the so-called compound network admittance 
matrix, (i.e., the 3-phase representation of the grid model which takes into account the various 
couplings between the network phases) instead of the single-phase equivalents. In this cases, 
each of the constraints in the OPF formulation needs to be formulated separately for each 
network phase. 
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A2. Any two-port component (e.g., lines, transformers etc.) is represented as 
a 7 r-equivalent; 

A3. We assume a perfect knowledge of the system parameters, i.e., the network 
admittance matrix is known; 

A4. The nodal-power injections are voltage-independent; 

A5. The control variables are composed by the nodal power injections/absorptions. 

2.1. The Proposed Centralized OPF Algorithm 

We are interested in maximizing the social welfare of the economic agents 
that use the grid, while maintaining an acceptable network voltage profile and 
respecting the line ampacity limits. Specifically, we tune the line ampacities and 
the network voltage profiles by controlling the (P, Q)-injections of distributed 
controllable devices Q (e.g., renewable generators) in a “fair” way: Each con¬ 
trollable device g £ Q has a certain utility function U g (•), and the sum of 
these utility functions is maximized subject to the satisfaction of the network 
operation constraints (voltage and ampacity). The resulting set-point is thus 
Pareto-optimal, i.e., no single device can increase its utility without hurting the 
utility of some other device, and locally-“fair”, i.e., the resulting set-point is a 
local maximizer of the sum of the device utilities lying on the Pareto boundary 
of feasible set-points. 

By convention, each line t £ C has a “receiving” and a “sending” end, which 
we denote by and t ~, respectively. These are chosen arbitrarily. A line is con¬ 
nected to two adjacent buses to which we refer by fd((. + ) and /3(£~), respectively. 

For each line, we introduce two auxiliary variables E^+ and E^~ representing 
the complex voltage at the two ends of the line. Assumptions A1-A3 allow us 
to express the corresponding injected currents and powers at the two ends of 
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line £: 


h+ — h+ (Ee+ 1 E e -) 

) = {Y t + Y ( +)E e + - Y ( E e - 

(1) 

h- = h- {Ee+, E e ~) 

1 = {Yt + Y e -)E e - — Y(E(+ 

(2) 

S e+ =S l+ {E l+ ,E e -) 

1 = Et+l e+ 

(3) 

S t - =S e -(E i+ ,E e -] 

1 = Ei- li~ 

(4) 


In the remainder of this paper, unless otherwise stated, the complex line currents 
and powers expressed above are always computed according to equations (l)-(4). 
They are thus all functions of Eg+ and E(- exclusively, although the arguments 
are often omitted for the sake of brevity. All quantities are expressed in “per- 
unit”, unless otherwise specified. 

For readability, we denote the vector formed by the real and imaginary 
parts of variables (Ei+,Ei-)i by y € R 4L , where L = \C\ is the number of lines. 
Note that for a given value of y, the corresponding currents and powers do not 
necessarily satisfy Kirchhoff’s law. 

We call y feasible if it satisfies voltage consistency and per-bus power- 
balance. Voltage consistency means that the voltages of all the lines incident 
to a specific bus b £ B are identical, i.e., have the same amplitude V& and the 
same argument ipi,: 

\E e+ \ = V 0{e+) , \E e -\=V 0{e - ) (5) 

arg {E e +) = tp&ig(E e -) = W € C. (6) 

At each bus b £ B, power-balance is satisfied if and only if 

]T -v+ ]T s e -=-Y,s g -m, v&gr, (7) 

P (<?+)=& 0{e~)=b g&b 

where S g is the controlled generated power of device g found at bus b, S(b) 
denotes the non-controllable power injection at bus b, and S (+, S(~ are obtained 
via (3)-(4). 

If y is feasible, it is important to note that equations (l)-(4) describe the 
exact AC power-flow equations. Hence, we use a non-approximated model of 
the grid. 
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We write the OPF formulation (Part I) equivalently: 2 


max \ Wg(Sg) subject to: 

s gi v b , Vb 

B B 9 66 

E e+ ’ E e- 


(8) 

Feasibility constraints (5), (6), (7) 



\I(+ | A, max aild \I e - | It'tnax, 

VI £C 

0) 

Vmin < 14 < Vinax, 

VbeB 

(10) 

Sg € Hg, 

v 5 ee 

(11) 


As previously stated, the objective function is the sum of the welfare of the 
controllable devices W g . In the above formulation, we denote by Q the set of 
controllable devices and by S g = P g + jQ q the controllable injected power by 
device g, subject to the capability constraint (11). The set Q can contain both 
generators and consumers. However, for the sake of presentation clarity, we con¬ 
sider that Q contains uniquely PV generators. This is not a limiting assumption, 
as our results apply to any device with controllable power injections (including 
controllable loads). Non-controllable loads do not appear in the objective func¬ 
tion, that expresses the utility of PV generators (a concave increasing function 
U(-) of active power injection) and the losses of the power converter: 

W g (S g ) = Ug(Pg) - T](Pg + Q 2 g), V<7 € 0. (12) 

We consider typical capability curves of PV power inverters: 

H g = {S g : \S g \ < ^,max5 \ ar 9{Sg)\ < <£ g, max }. (13) 

In order to solve the problem (8)-(ll), we convert the inequality constraints 
(9) to equality constraints by introducing slack variables if+ and ig- as follows: 

It?+1 + H+ — h ,max and \Ig-\ + i t - = Ie ,max) W G C (14) 

i e +,it-> 0, (15) 


2 Unlike in Part I, we consider wlog that there are two types of connected devices: they 
either have controllable power injection S g or impose an overall fixed power injection S(b) in 
bus b. 
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We denote by x the real vector of variables formed by the artificial control 
variables (14, (fib)beB, {H+Ae.-)leCi and the device controllable injected power 
(Pgi Qg)geS- 

Notice that all the equality constraints above, (5), (6), (7), and (14) can 
be summarized as g(y) + Ax + b = 0, where g(-) is a smooth non-convex func¬ 
tion that can be derived from equations (l)-(6), and A is a positive definite 
matrix. Similarly, the inequality constraints, (10), (11), and (15), can be ex¬ 
pressed as h(x) > 0, where h(x) is a convex function that can be derived from 
equations (10), (15), and (13). We denote the objective by f(x), where / is 
concave. 

We can thus write our problem in the more compact form: 


maxf(x) 

(16) 

x,y 

subject to g(y) + Ax + 6 = 0 

(17) 

h(x) > 0. 

(18) 


We write its augmented Lagrangian ([9, 10, 11]): 

£p{x, y ; A) =f(x) + A !{g{y) + Ax + b) 

- P -\\g(y)+Ax + b\W (19) 

where p is the weight of the quadratic penalty term added to the classic La¬ 
grangian function, and A is the vector of Lagrange multipliers associated with 
the equality constraints (17). 

Our centralized iterative algorithm for solving the OPF is based on the 
method of multipliers ([9, §4.2]). This method was first introduced for solving 
iteratively non-linear equality constrained problems. It is shown to converge un¬ 
der more general conditions than dual ascent [22]. Algorithm 1 summarizes the 
proposed centralized algorithm, and Theorem 1 characterizes its convergence. 

The main advantage of the method of multipliers is that there exists a finite 
value p such that the problem (20) is locally convex for all p k > p. Note also 
that the algorithm bounds the value of A at each iteration. The next vector of 



multiplier estimates A is obtained after a projection on the set [—A, A] defined as 
[— Ai, Ai] x [—A 2 , A 2 ] x ...; the constant vector A is chosen such that the sought 
optimal vector of Lagrange multipliers A* lies in [—A, A] (see [23, §2.2.2]). 

Theorem 1. For smooth objective function f £ C 2 and suitably chosen A such 
that the optimal vector of Lagrange multipliers X* satisfies X* £ [—A, A], Algo¬ 
rithm 1 converges to a local minimum of the nonlinear progam (16)-(18). 

Proof. By [23, Proposition 1.23], our problem satisfies assumption (S) from [23, 
§2.2], since the equality constraint is a C 2 function of y , and the objective func¬ 
tion is chosen to be C 2 . Proposition 2.7 from the same reference guarantees the 
desired convergence, if the iterates (x k ,y k , X k ) reach the set D from Proposition 
2.4 of [23], i.e., if there exists a k such that ( x k , y k , X k ) £ D (for all the following 


Algorithm 1 Centralized algorithm for the OPF (16)-(18) 

• Set k= 0 and initialize control variables x and y: 

= 0,£$+ = E\_ = l,U b ° = l,<pg = 0 ,i\+ = i°_ = 0 (per-unit), 
Lagrange multipliers A 0 = 0, increasing gain sequence ( p k )k , p k — > 00 . 

1: repeat 

2: Maximize the augmented Lagrangian for fixed A = X k : 

(a; fc+1 , 2 / fc+1 ) = arg max L pk {x,y\X k ). (20) 

x,y:h(x)> 0 

3: Update the Lagrange multipliers: 

A fc+1 = {A fc + p k [g(y k+1 ) + Ax k+1 + b] } (21) 

4: k <r- k+1 

5: until the maximum number of iterations has been reached or the change 
in the Lagrange multipliers between two consecutive iterations is less than 
a tolerance S > 0 
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indices k > k, the iterates stay in D , and convergence ensues). The existence 
of such a k follows from the choice of the divergent increasing sequence of gains 
( p k ) and from the boundedness of the sequence (A fc ). □ 

Due to the quadratic terms in the expression of the augmented Lagrangian (19), 
the optimization problem in (20) does not decouple across the network and, 
therefore, cannot be solved in a distributed manner. In the following section, 
we reformulate this problem in an equivalent way that leads to a distributed 
algorithm for its resolution. 

2.2. Distributed Solution of the OPF Problem 

We adopt a primal decomposition method [21] that gives an iterative algo¬ 
rithm for the minimization of the problem in Step 2 of Algorithm 1. In (19) 
the line voltages y = (£^+, E^~ ) are “coupling” variables. If these variables are 
fixed to a specific value, then problem (20) decouples in smaller, easier (convex) 
problems, that can be solved by local agents. 

Specifically, to solve (20) iteratively for fixed values of the Lagrange multi¬ 
plier estimates A and fixed gain p we take the following approach: At the n-th 
iteration, the value of the coupling variables y n = (E™ + ,E™_) is assumed fixed. 
The x variables, i.e., the power set-points of the controllable devices ( S g ), the 
bus voltages (Vb), and the slack variables i(+,i^~, are computed by solving the 
following constrained convex optimization problem: 


x n+1 = arg maxLp(j), y n , A). (22) 

x:h(x)> 0 

Next, the coupling variables y are updated as follows: 

V n+1 =V n + a n (V y L f> )(x n+1 ,y n , A), (23) 

where a n is a positive step-size sequence of the gradient descent. The choice 
of the step-size is related to the topology of the network and the parameters of 
the lines (i.e., the network admittance matrix). For example, a large constant 


10 



step-size might not allow the algorithm to converge, whereas a small constant 
step-size could cause slow convergence 3 . 

The algorithm stops when the norm of the update in the y variables is less 
than some small positive tolerance e, i.e when || X/ y Lp(x n+1 , y n , A)|| < e. 

Theorem 2. The algorithm (22)-(23) with tolerance e in the stopping criterion 
converges to a vicinity B{(x* ,y*), S) of a local optimum (x*,y*) of problem (20). 
If (20) is strongly locally convex in y in a vicinity of (x*,y*), then 6 = 0(e 2 ). 

Proof. (Sketch) Denote v(y) = max^.^^pQ Lp(x, y, A) and x*(y) the value of x 
that achieves this maximum (22). Theorem 2.1 of [24] says that the optimum 
(x*(y*),y*) of max y v(y) coincides with the one of (20). Moreover, a (5-optimal 
solution (x*(y$),ys) of ma x y v(y) (that is, v(ys) > v{y*) — S ) is also (5-optimal 
for (20). 

We now show that it holds that V y v(y) = (y y Lp)(x*(y),y, A), or equiva¬ 
lently, Dx d ^ j ' > (' V x Lp)(x*(y ), y, A) = 0. If we can show this, then the algorithm 
(22)-(23) is equivalent to a gradient ascent in y on v(y). It is easy to show 
that the function v(y) is “smooth” (C 2 ). By the strong local convexity around 
(x*,y*) of the augmented Lagrangian, [9, Exercise 1.2.10] allows us to conclude 
that S = 0(s 2 ). 

Note that problem (22) is convex. Consider the optimal multipliers y* cor¬ 
responding to the constraints h(x) > 0. They satisfy the KKT conditions: 

(V x Lp)(x*(y),y, A) = V*{y)V x h i( x * (y)) 

i 

T*(y) h i( x *(y)) = °; m *> o . 

Define the following functions: i)i{y) := hi(x*(y)). Since x*(y) is always feasible, 
it means that 4>i(y) > 0. Consider the set of indices To (y) ■= {* : hi(x*{y)) = 0}. 
Take some i € In this case the function ipi(y) l ias an extremal point in 

y, which implies that \/ y ifi(y) = 0, or again that — X D ^ S7 x hi(x*(y)) = 0. For 

3 In order to properly tune this parameter, a dedicated off-line study can be performed 
before deployment of the proposed algorithm. 
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all i ^ I\ } (y ), by KKT we have y*\y) = 0. By the above arguments, 

Dr* f ^ 

— fr^-(^xLp)(x*(y),y,X) 

Dy 

= X] (y) Vx/N (a:* (y)) = 0. 

□ 

Thanks to its separability property, problem (22) can be solved in a dis¬ 
tributed manner. Bus agents can be responsible for updating the power set- 
points of the controllable devices (S g ) that are connected to them, as well as 
their voltages (14) in parallel, and lines can be responsible for updating the 
slack variables (i(+ ,i{-)- Specifically, the ‘power set-points (<§” +1 ) of devices in 
bus b are obtained by solving the following convex problem: 

(5^ +1 ) = arg max^ W g (S g ) 

9&b /3(t+)=6 0(t~)=b ^ 

where Ah is the given multiplier corresponding to the constraint (7) of bus b. 
The other problems (for the other x variables) have simpler expressions that we 
do not reproduce for brevity sake. 

Similarly, (23) can be decomposed across the different network lines: line- 
agents can update the voltages at their two ends in parallel. In terms of required 
information, each bus agent needs to know only the voltage values of the lines 
that are incident to it, the constraints of the devices, and the state of the loads 
that are connected to it. Finally, in order to compute the partial derivatives 
of (23) with respect to its voltages, each line requires solely the information of 
the power balance and the voltage values of its two adjacent buses. The actual 
implementation of the distributed synchronous OPF algorithm is summarized 
below in Algorithm 2. 

Theorem 3. For smooth objective function f £ C 2 and suitably chosen X such 
that the optimal vector of Lagrange multipliers A* satisfies X* £ [—A, A], Algo¬ 
rithm 2 converges to a local minimum of the nonlinear progam (16). 
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Proof. The proof is similar to the one of Theorem 1. It uses Proposition 2.16 of 
[23] for convergence, which only requires at each iteration a (5 fe -optimal solution 
for (20) with 5 k —> 0. By Theorem 2 we can conclude. □ 

In a realistic setting, in order to take full advantage of the distributed for¬ 
mulation of the OPF algorithm, as described above and to avoid the overhead 
cost of coordination between agents, the updates should be performed in an 
asynchronous fashion. Contrary to ADMM-based algorithms, which require a 
synchronized implementation of the updates, the proposed algorithm can be im¬ 
plemented in an asynchronous manner. In this direction, we assume that each 
of the bus and line agents has its own two local poisson clocks with different 


Algorithm 2 Distributed algorithm for the OPF (16)-(18) 

• Set k= 0 and initialize control variables x and y: 

S° = 0 ,E° + = E° t _ = l,P b 0 = l,vl = 0,i° + = i° e _ = 0 (per-unit), 
Lagrange multipliers A 0 = 0, increasing diverging gain sequence ( p k )k , 
p k — > oo, decreasing tolerance sequence ( e k > 0)*, £ k —> 0. 

1: repeat 

2 : n «— 0; x° «— x k ; y° <— y k 

3: repeat 

4: x n+1 = arg max x . h{x) > 0 L pk (. x , y n , X k ) 

5: y n+1 = y n + a n (V y L pk )(x n+1 ,y n ,\ k ) 

6: Tl <r- n+l 

7: until \\S7yL p k (x n+1 ,y n , A fc )|| < e k 

8: x k+1 t- x n+l ; y k+1 <- y n+1 

9: A fe+1 = n [ _ X S] {A fc + P k [g(y k+1 ) + Ax k+1 + b] } 

10: k t- k +1 

11: until the maximum number of iterations has been reached or the change 
in the Lagrange multipliers between two consecutive iterations is less than 
a tolerance S > 0 
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rates. The clock with the lower rate (C'l) triggers the multiplier update (21) 
and the clock with the higher rate (C 2) triggers the events described in steps 
(22)-(23). 

In detail, all the control variables and the Lagrange multipliers are first 
initialized. Then, every time the C 2 clock of a bus ticks, this bus performs 
local update operations by using the most recent stored values for the voltages 
of its incident lines and for the associated Lagrange multipliers. Once the bus 
updates its power and voltage values, it informs the incident lines of the changes. 
Similarly, when the C 2 clock of a line ticks, the line agent updates the variables 
( i /+, i(-) by taking into account the most recent values of the line current flows 
and associated Lagrange multipliers. In addition to this update, the updates of 
the voltages of its two end-points are triggered. In order to compute the new 
values, the line uses the most recent stored values for the adjacent buses’ powers 
and voltages, and once the updates are completed the line communicates this 
information to its neighboring buses. Now, when the C'l clock of a bus or a line 
ticks, then the corresponding agent updates the Lagrange multipliers (21). It 
is worth noting, that we no longer have a serial implementation of the various 
updates like the ones presented in Algorithm 2. On the contrary, the different 
rates of the clocks are chosen in such a way to ensure that, on average, a sufficient 
number of the updates occurs before an update of the corresponding Lagrange 
multiplier takes place. 

In this section, we investigate the performances and convergence properties 
of the centralized Algorithm 1 in several different scenarios. In particular, we 
consider the cases presented in Part I of the paper, where the BFM convex- 
ification leads to an incorrect solution of the OPF problem and ADMM fails 
to converge to a solution. Additionally, we investigate the performances of the 
proposed centralized algorithm under different initial conditions of the electrical- 
network state. In order to do so, we consider the same 4-bus test network that 
was used in Part I of the paper. We assume a first test case where the con¬ 
trollable device connected to bus 4 is a generator, whereas controllable loads 
are connected to buses 2 and 3. The network characteristics, the base values, 
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the capability limits of the controllable resources, and the voltage and ampacity 
limits are given in Part I (Fig. 3 and Table 7). In what follows, the objective 
function accounts for the minimization of the network losses, as well as for the 
utility of the generating units, namely: 

min - 'S'' Re(Sg) + Re(Y e )\R \ 2 (24) 

2.3. Effect of the Line Length, Network Rated Value and Network State on the 
Convergence of Algorithm, 1 

In order to compare the performances of the proposed algorithm with the 
OPF algorithm proposed in [1, 2], we solve the OPF problem for various line 
lengths and network voltage rated values as in Part I of the paper. In particular, 
we assume that the line lengths are uniformly multiplied by a factor in the range 
[1.25 — 7.5] (while keeping the network voltage rated value to its nominal value) 
and the network voltage rated value varies in the range [15 — 40 \kV (while 
keeping the line lengths to their nominal values). The evolution of the bus 
voltages, the line-current flows, as well as the active and reactive powers, are 
shown in Figures 1-6. It is worth noting that in all the cases the proposed 
algorithm converges in a few iterations. Furthermore, we observe from Fig. 2 
and Fig. 5 that the line-current flows satisfy the line ampacity limit, once the 
algorithm has converged, in all cases. In particular, in Fig. 2 it is worth observing 
that as the line length increases the receiving and sending-end current flows of 
the same line become significantly different. The behavior of the current flows 
as the voltage rated value increases is similar (Fig. 5). This effect is due to 
the increasing contribution of the current flow toward the shunt elements of 
the lines. In fact, we show, in Figures 7 and 8, the amount of reactive power 
produced by the shunt elements of the lines for the various values of the line 
lengths and the network voltage rated values. We observe that as the line length 
increases or the rated value of the voltage increases the reactive power produced 
by the shunt elements of the line increases as well. 

We investigate, in addition to the effect of the line lengths and the network 
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Fig. 1: Evolution of the magnitude of network voltages for various line lengths. 



Fig. 2: Evolution of the line current flows for various line lengths. 


voltage rated value, the performance of the proposed algorithm under a different 
network operating point. To this end, we consider a second test-case where the 
controllable device connected to bus 4 is a controllable load and generators are 
connected to buses 2 and 3. In this respect, we consider an extra term in the 
objective function, which represents the utility associated with the controllable 
load and is given by (Pl — P 0 ) 2 , where P 0 represents a constant amount of load 
that has to be served. The capability limits of the controllable resources are 
shown in Table 1. The convergence of the voltages, current flows, as well as 
active and reactive powers are shown in Fig. 9. For the sake of brevity, we only 
show the evolution of the active and reactive power of the controllable load of 
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Fig. 3: Evolution of the active and reactive power of the controllable devices for various line 
lengths. 
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Fig. 4: Evolution of the magnitude of network voltages for various values of the network rated 
voltage. 


bus 4, as the controllable generators are small and reach their maximum value 
upon convergence. 

2.f- Performance Evaluation of the Proposed Algorithm in the Presence of Shunt 
Capacitors in the Network 

In what follows, we consider the same network adopted in the previous sec¬ 
tion and a case where each network bus, apart from the slack, has a load and 
a generator connected to it. The demand in the network is assumed to be non- 
controllable, whereas the generators are assumed to be distributed solar panels 
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Fig. 5: Evolution of the line current flows for various values of the network rated voltage. 



Fig. 6: Evolution of the magnitude of active and reactive power of the controllable devices for 
various values of the network rated voltage. 
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Fig. 7: Reactive power produced by the shunt elements of the lines for various values of the 
network voltage rated value. 
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Fig. 8: Reactive power produced by the shunt elements of the lines for various line lengths. 


Table 1: Parameters of the test network used for the investigation of the performance of the 
proposed OPF algorithm under a different operating point 


Parameter 

value 

[P 9 min ,P gm J(bus 2) (MW) 

[0,0.01] 

[P Bmtn ,P Bm J(bus 3) (MW) 

[0,0.012] 

(Pc min ,Qc min )(M W,Mvar) (bus 4) 

0.3,0.15 

P 0 (MW) (bus 4) 

1 



Fig. 9: Evolution of the magnitude of network voltages, current flows, as well as active and 
reactive power of the controllable load at bus 4 for the case of low generation and high load 
in the network. 
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Table 2: Parameters of the test network used for the evaluation of Algorithm 1 in the presence 
of shunt capacitors in the network 


Parameter 

Value 

Generators’ power, \Si gmax A = 2, 3,4 (MVA) 

0.40,0.39,0.46 

Generators’ power factor, cos<pi g ,i = 2,3,4 

0.9 

Loads’ active power, P, c ,* = 2,3,4 (MW) 

2.76,2.16,2.46 

Loads’ reactive power, Qi c ,i = 2,3,4 (MW) 

1.38,1.08,1.23 

Shunt capacitor (bus 2)(uF) 

859 

Penalty term gain, p 

10 4 

Tolerance and maximum number of iterations 

KT 4 ,10 4 

[^ram? ^matc] (p*u) 

[0.9,1.1] 


with typical PV-type capability constraints given by (13). For this scenario, the 
capability limits and the values of loads and generation are shown in Table 2. 
In addition to the loads and generation, we consider that a shunt capacitor is 
connected to bus 2. In order to model this shunt capacitor, we consider that it 
is part of the first line. In particular, we consider that the shunt capacitance 
on the sending end of the 7r-model of the line that connects buses 1 and 2 is 
modified accordingly, to account for the shunt capacitor. For this particular test 
case, it is worth noting that ADMM exhibits oscillations and fails to converge 
to a solution (see Part I, Fig.9-12). 

The results for this specific test-case, for the voltage magnitudes and the 
active and reactive power of the buses, are shown in Fig. 10. It is worth ob¬ 
serving that the proposed algorithm converges to a solution within a few tens of 
iterations; which is contrary to the ADMM-based solution of the OPF problem. 

2.5. Coordinated Control of DERs Power Set-points and OLTC 

In this section we investigate the performances of the centralized OPF al¬ 
gorithm in the case of coordinated control of the DERs’ nodal power injections 
and the OLTCs positions. We assume that at the primary substation there 
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Number of iterations Number of iterations 


Fig. 10: Evolution of the active and reactive power, as well as the voltages of the buses when 
a shunt capacitor is connected to bus 2. 


is an OLTC that can vary the secondary voltage within a ±6% range in ±36 
steps. In order to do so, we assume as described in Section 4.2 of Part I of the 
paper that the first network branch connecting network nodes 1 and 2 has a 
longitudinal admittance that is a function of the OLTC ratio a and the trans¬ 
former nominal admittance, namely aY t . The shunt elements of the receiving 
and sending ends of this specific branch are also a function of the new control 
variable, i.e., a(a — \)Y t and (1 — a)Y t respectively (based on the OLTC model 
in [25]). We model this control variable a £ I as pseudo-continuous and round 
it to the nearest integer when the algorithm has converged to a solution. We 
also include a cost function in the objective function of the problem related to 
the OLTC operation. In particular, the changes of OLTC need to be penalized 
as these devices are typically used by the DNO rarely due to their increased cost 
and their limited lifetime. This is why we include a quadratic penalty function 
at each time-step in order to minimize the OLTC changes from their previous 
position. In the application example that follows, we consider the initial OLTC 
position set to -18, the rest of the grid parameters are as in Table 2 except for 
the voltage limits that are set in this case to ±5% of the network rated value. 
The results for this test case are shown in Fig. 11 and 12 below. In particular, 
Fig. 11 shows the evolution of the voltage magnitudes in all network buses and 
in solid black line the evolution of the control variable a, which upon solution of 
the problem is rounded to position -4. Fig. 12 shows the evolution of the rest of 
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Fig. 11: Evolution of the magnitude of network voltages when OLTC is included as a control 
variable. 



Fig. 12: Evolution of the magnitude of the nodal active and reactive power of the network 
buses when OLTC is included as a control variable. 


the control variables, which are essentially the nodal active and reactive power 
injections. Overall we can observe that for this test case the OPF algorithm 
converges in 31 iterations and the OLTC ratio converges to the value which is 
closest to its initial setting and at the same time allows the voltage magnitudes 
at all buses to lie within the acceptable limits for safe operation. 
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Fig. 13: Evolution of the magnitude of network voltages, line current flows and active and 
reactive power of the controllable devices when the initial voltage magnitudes are set to 0.9 
and the voltage angles to —7r/6. 


2.6. Performance Evaluation of the Proposed Algorithm under Different Initial 
Conditions of the Network State 

Finally, we investigate the performances of the proposed algorithm under 
different initial conditions of the network state variables. In order to do so, we 
initialize the magnitude of the control variables E ® + , Eff _, V)° in Algorithm 1 in 
the range [0.9,1.1] and their angle in the range [—7r/6, 7 t/ 6], totaling 121 different 
cases. For each combination, we solve the centralized OPF problem for the same 
network adopted in Part I (Fig. 3). In all the cases the algorithm converges to 
the same solution within a few tens of iterations. In Table 3, the mean value of 
the number of iterations, as well as the 95-th percentile are shown. For the sake 
of brevity, we show in Fig. 13-14 the convergence results for the voltage, as well 
as for the current flows and the active and reactive power profiles for the two 
extreme cases, specifically when the voltage magnitude is set to 0.9 (1.1) and 
the voltage angle is set to —7 t/6 (7t/6). 


Table 3: Number of iterations for the solution of the OPF problem (Algorithm 1) 



Mean number of iterations 

95-th Percentile 

Algorithm 1 

18.21 

46.45 
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Fig. 14: Evolution of the magnitude of network voltages, line current flows and active and 
reactive power of the controllable devices when the initial voltage magnitudes are set to 1.1 
and the voltage angles to 7 t/6. 


3. Performance Evaluation of the Proposed Distributed Asynchronous 

OPF Algorithm 

For the sake of completeness, in this section, we assess the performance of 
the proposed algorithm with respect to a realistic grid represented by a modified 
IEEE 13-node test feeder ([26]). The modifications are (i) balanced lines, (ii) 
inclusion of secondary substations where voltage independent PQ-injections are 
placed, and (iii) lines ten times longer. We use this benchmark to assess the 
behavior of the proposed distributed asynchronous OPF algorithm. Also, we 
compare the solution and convergence of the distributed version of the algorithm 
to the centralized one. 

We consider a test case where each network bus, apart from the slack bus, has 
a load and a generator connected to it. The demand in the network is assumed 
to be non-controllable, whereas the generators are assumed to be distributed 
solar panels with typical PV-type capability constraints. For this test case, the 
capability limits and the values of loads and generation are shown in Table 4. 

We solve the OPF problem in (8)-(13) using Algorithm 1, as well as the 
asynchronous implementation of Algorithm 2. The results are shown in Fig. 15- 
17. For the sake of brevity, we plot only the evolution of the magnitudes of the 
minimum voltage, the maximum voltage and the median value of the voltage. 
We plot also the evolution of the minimum, maximum and mean values of the 
current flows on the receiving-end of the line and the evolution of the active and 
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Table 4: Capability limits and values of loads and generation for the evaluation of Algorithm 
2 _ 


Bus 

s gmax 

(MVA) 

P C (MW)/ 

Q c ( Mvar) 

Bus 

s 9max 

(MVA) 

P c { MW)/ 

Q c (Mvar) 

2 

0.0437 

0.0025 / 0.0011 

8 

0.0347 

0.0031 / 0.0014 

3 

0.0480 

0.0029 / 0.0012 

9 

0.0403 

0.0031 / 0.0013 

4 

0.0506 

0.0032 / 0.0013 

10 

0.0373 

0.0031 / 0.0013 

5 

0.0367 

0.0029 / 0.0012 

11 

0.0482 

0.0024 / 0.0010 

6 

0.0443 

0.0029 / 0.0012 

12 

0.0399 

0.0030 / 0.0013 

7 

0.0426 

0.0025 / 0.0010 

13 

0.0436 

0.0029 / 0.0012 


reactive powers. It is worth observing that Algorithm 1 converges to the opti¬ 
mal solution within a few iterations and also that the distributed asynchronous 
implementation of Algorithm 1 converges to the same solution as its centralized 
counterpart. 




Fig. 15: Evolution of the voltage magnitude for the distributed asynchronous algorithm as a 
function of the number of messages exchanged (left) and for Algorithm 1 as a function of the 
number of iterations (right). 
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Fig. 16: Evolution of the current flows for the distributed asynchronous algorithm as a function 
of the number of messages exchanged (left) and for Algorithm 1 as a function of the number 
of iterations (right). 
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Fig. 17: Evolution of the active and reactive power for the distributed asynchronous algorithm 
as a function of the number of messages exchanged (left) and for Algorithm 1 as a function 
of the number of iterations (right). 


4. Conclusion 

To overcome the limitations identified in Part I, we have proposed algorithms 
for the solution of the AC non-convex OPF problem in radial networks that are 
proven to converge to a local minimum. These algorithms use an augmented 
Lagrangian approach and rely on the method of multipliers for the OPF solu¬ 
tion. The two algorithms solve the centralized and decentralized (asynchronous) 
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formulation of the targeted OPF. We have shown the robustness of the central¬ 
ized version with respect to the following elements: (i) various line lengths, (ii) 
various network-rated voltage values and (ii) different network operating points 
(cases where the BFM convexification leads to an incorrect solution), (iii) the 
presence of shunt capacitors in the grid (where ADMM failed to converge to 
a solution) and (iv) different initial conditions of the electrical network state. 
Finally, we have verified the equivalence of the two proposed algorithms for 
the case of the IEEE 13-node test distribution feeder where realistic operating 
conditions have been considered. 
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